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We study the backreaction of free quantum fields on a flat Robertson- Walker spacetime. Apart 
from renormalization freedom, the vacuum energy receives contributions from both the trace 
anomaly and the thermal nature of the quantum state. The former represents a dynamical realisa- 
tion of dark energy, while the latter mimics an effective dark matter component. The semiclassical 
dynamics yield two classes of asymptotically stable solutions. The first reproduces the ACDM model 
in a suitable regime. The second lacks a classical counterpart, but is in excellent agreement with 
recent observations. 
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Introduction. During the last decade, cosmological obser- 
vations have established an ongoing phase of accelerated 
expansion in the recent history of our Universe [l|-l5| . Ac- 
cording to the ACDM concordance model of cosmology, 
the acceleration results from a modification of Einstein's 
equations by a cosmological constant (A) which has just 
begun to dominate the total energy density of the Uni- 
verse. More generally, such a contribution is called dark 
energy. In addition, the Universe is widely believed to 
be filled with cold dark matter (CDM), a non-luminous 
and weakly interacting component responsible for the for- 
mation and growth of large scale structures [6|48| , which 
outweighs the contribution of baryonic matter by a fac- 
tor of ten. Despite all efforts, the dark energy lacks a 
sound theoretical understanding so far. The purpose of 
this letter is to shed a new light on these issues from the 
fundamental point of view of quantum field theory on 
curved spacetimes. 

Compared to earlier attempts on this task - see for 
example [9l4l3| - the results presented here benefit from 
recent progress on a quantization scheme which is intrin- 
sically tied to non-trivial spacetimes [lj, LL5J , and repre- 
sents a significant leap forward in our understanding of 
quantum field theory on curved backgrounds. Our aim is 
to apply these new techniques in a cosmological frame- 
work in order to address the issue whether dark energy 
and, if possible, dark matter originate from fundamental, 
non-interacting quantum fields. 

Cosmological quantum fields. We consider a four- 
dimensional flat Robertson- Walker (RW) spacetime, 
whose metric ds 2 = —di 2 + a 2 (t)dx is completely deter- 
mined by the scale factor a(t). The dynamical evolution 
of the Universe in terms of a(t) is to be obtained from a 
solution of the semiclassical Einstein equations, 

G^v — 87rG(:T /iI/ :)n , 

where G M „ is the Einstein tensor and G the gravitational 
constant. On the right hand side, one has to take the 
expectation value of the regularised matter stress-energy 



tensor in a Gaussian state f2. The latter is chosen to ful- 
fil the so-called Hadamard condition which fixes the UV 
properties of Q and assures the existence of a well-defined 
notion of normal ordering of fields, : • : [lj]. The consis- 
tency requirement that the expectation value of {'Ty, v :)n 
should be covariantly conserved unavoidably leads to a 
term in the trace of (:T M „ :)o which has no classical coun- 
terpart, the trace anomaly [16Hl8( |. Given a collection of 
A^o conformally coupled real scalar fields <j>i, i = 1, ...,No, 
with mass rm, iVi/2 Dirac spinors ipj, j = 1, ..., Ni/ 2 with 
mass rrij, and N\ massless vector fields on a RW space- 



time, the trace reads 19142 



(: T^ :) n = a (rTR^ - ^- j + a{m)R + e 2 (m) + 

N Wl/2 

+c 3 DR - J2 m U-ti -)n - Y. m i^^i : >n ■ (1) 

Here, R^ is the Ricci tensor and R the scalar curvature, 
while a = {2S%Qii 2 )- 1 {N + llN 1/2 + 62Ni). Moreover, 
c\ (to) and c 2 (to) are linear combinations of powers of the 
masses. The corresponding mass coefficients, as well as 
the coefficient C3 itself, are subject to a finite renormal- 
ization freedom, and hence free parameters of the theory. 
Note that any state-dependent quantity appears multi- 
plied by a suitable power of the mass, so that massless 
fields contribute only to the purely geometric terms. 

Observing the homogeneity and isotropy of our met- 
ric ansatz, the semiclassical Einstein equations can be 
rewritten in terms of two coupled inhomogeneous ordi- 
nary differential equations, 
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where dots denote derivatives with respect to t. Here, 
H = ^ is the Hubble function and we set p = (:Too -)n- 
The quantum energy density p is obtained from ^j up 
to a solution of the corresponding homogeneous equation, 
which is of the form p = Coa -4 . As we are ultimately 



interested in describing the late time evolution of the Uni- 
verse, which is characterised by a low spacetime curva- 
ture, we shall employ an "adiabatic" approximation and 
discard terms in p with time derivatives of H. Hence, the 
contribution of the purely geometric, state-independent 
part in ((T|) to the quantum energy density reads 



Pgoom = 3aH 4 - 3ci(m)H 2 - 



c 2 {m) 



As H <C m for generic masses in the late cosmic history, 
we supplement the adiabatic approximation by discard- 
ing 0(H 5 ) terms in p. We retain 0(H A ) terms to cap- 
ture the lowest order effects of curved spacetime quantum 
fields. 

State dependence. We consider every quantum field to be 
in a state VL which fulfils an approximate KMS condition 
at some instant of time t = to in the past [23]. Hence, 
fi depends on a fixed temperature parameter T equal for 
every field. Moreover, we demand that T <^ m and com- 
pute all terms to the lowest non-trivial order in T/m. We 
start with the conformally coupled Klein-Gordon field. 
Following 23j, \2A | , the two-point function of Q reads 
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where k = \k\ and 4>t{%) = Tfc(t)e are suitably nor- 
malized Klein-Gordon modes. Furthermore, fco = {k 2 + 
rn2fl o) 1 an( i a o = a(to), whereas /3 is the inverse tem- 
perature. On account of our choice of : • : and $7, and up 
to the renormalization freedom already present in (JTJ , we 
find 
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27r 2 [ e /3fe _ 1 l Tfe WI ^ k 2 + m 2 a 2 
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Considering our adiabatic approximation, we compute 
|7fe(t)| 2 by expanding it in terms of H and its time deriva- 
tives so that the resulting p depends only on H n , n < 5. 
This yields 

\T k (t)\ 2 =91 + 92H 2 + g 3 H + g A H 4 + g 5 HH 2 . 

Here, gi for i = 1, ..., 5 are suitable functions of a, k, and 
m. Since |7\.(rj)| 2 fulfils an ordinary differential equation 
which descends from the Klein-Gordon equation, these 
functions can be determined recursively by taking the 
UV properties of O as an initial datum. As shown in 
detail in 



241 ] . we thus obtain 
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up to the renormalization freedom accounted for in ([T]) 
and up to higher orders in T/m. Here, / is a monotoni- 
cally decreasing function of aofim with 1(0) = C(3). 



For Dirac fields, the object of interest is the expecta- 
tion value (:tp^p:)n- As explained in 23|,|24|, and barring 
the renormalization freedom already present in ([T]), this 
quantity can be computed as 
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where Qfc are normalised Dirac modes. Recalling our 
adiabatic approximation, we expand |Qfe| 2 in terms of H 
and its derivatives up to sufficient powers. This yields 

\Q k \ 2 = Gi + G 2 H + G 3 H 2 + dHH 2 + G 5 H 4 , 

where again the Gi depend on a, m and k. By a proce- 
dure similar to the scalar case we obtain 
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19H 4 - 15Hm 2 + 17 HH 2 
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see [24j for the details. Here, L is a function sharing the 
properties of /. 

Collecting the state dependent terms of (:T M „ :)n and 
computing the corresponding energy density, one finds 
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Here, 04(771) > is a fixed and constant linear combina- 
tion of the field masses, and derivatives of H have been 
discarded in accord with our adiabatic approximation. 

Effective Friedmann equation. We insert p — p seo m + 
Pstatc into the second equation in ([2]) to obtain an alge- 
braic equation for H . If N x / 2 > 0, this equation is of 
fourth order in H. Its solutions may be interpreted as 
effective Friedmann equations and read 



H\a)± = Hi±\H$ 
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where if* depends on the number of fields and the renor- 
malization freedom, C 2 depends on the field masses and 
T, C3 is a free renormalization parameter, and G\ is a 
multiple of the already mentioned integration constant 
Co present in p. A qualitative analysis of ([3]) displays 
two asymptotically stable fixed points, corresponding to 
de Sitter spacetime. Physical solutions approach one of 
these points with H 2 > and H < 0, see figure Q] The 
value of the effective cosmological constant associated to 
the attractors is not fixed since it depends on the renor- 
malization constants of the underlying quantum field the- 
ory E|. 



Upper branch. Considering H+(a), the effective cosmo- 
logical constant at the fixed point turns out to be de- 
termined by the free parameters H* and C3, whereas 




FIG. 1: Generic plot of the two branches. While H+ displays 
a Big Bang singularity, H- becomes imaginary at finite time 
[22| . The latter phenomenon is an artefact of the adiabatic 
approximation employed. 
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FIG. 2: Best fits of the type la supernova distance moduli /i 
in the Union2 compilation [j. The solid line corresponds to 
the best fit upper branch, the dashed line corresponds to the 
best fit ACDM model or the lower branch. 
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FIG. 3: The Om(z) diagnostic of the upper branch (solid line) 
in comparison with the one of the ACDM model or the lower 
branch (dashed line). 



the physical interpretation of the remaining terms is far 
from obvious. In particular, there is no contribution from 
quantum fields behaving like (dark) matter or radiation. 
However, the upper branch solution describes the late- 
time evolution of the Universe at least as well as the 
concordance model does. In order to demonstrate this, 
we consider H + in terms of the redshift z = a^ 1 — 1 and 
perform a x 2 ~fit °f the latest type la supernova data, 
the Union2 compilation [4] , cf. figure O To empha- 
sise the fundamental deviation of the upper branch from 
the ACDM model, we set C3 = by hand and find the 
best-fit x 2 to be as low as the one of the ACDM model, 
namely, Xmin/^ / ~ 0.98. The best fit parameters are 
H* = 0.57flffo , Ci = -0.273^ , C 2 = 0.270F^, 
where H = H + (z — 0). The discrepancy from the 
predictions of ACDM can be made evident if we com- 
pare the two different models by means of the Om(z) = 
\(H(z)/Hn) 2 - 1]/[(1 + z) 3 - 1] diagnostic recently intro- 
duced in [25], see figure [3] This observable is particularly 
well-suited to determine whether or not dark energy is a 
cosmological constant at low redshifts. 

To summarise, the cosmological solutions associated 
with the upper branch provide a genuinely quantum the- 
oretical realisation of dynamical dark energy which can 
easily be distinguished from a cosmological constant by 
upcoming cosmological observations. However, it re- 
mains unclear how to reconcile the cosmological evolu- 
tion with the standard picture of matter and radiation 
domination in earlier stages of cosmic history. 

Lower branch. Let us now consider the class of solutions 
dubbed H- (a) in ([3|) . For sufficiently large values of H* , 
we can expand the square root to approximate H_(a) as 
Ki+K2d~ 3 +K3a~ 6 , where the constants Ki are specified 
in terms of C2, C3 and H* and where we have set C\ = 0. 
Since the value of K3 turns out to be of order 10~ 3 Hq 
in the best fit result, the model provides an effective de- 
scription of the concordance model, reproducing all the 
observational features of ACDM in the recent past. How- 
ever, deviations emerging from the higher order correc- 
tions are to be expected in the z>l regime, which is not 
yet observationally accessible. Moreover, the fit implies 

that C2, which equals (X^=°i rn i + 3j2jl{ 2 m j)T 3 U P to 
O(l) numerical factors, must be of the same order as the 
critical density ZH^^tiG)^ 1 . The fit is not sufficient to 
fix the masses and T at the same time. However, if we 
assume a single mass parameter m for simplicity, we get 
TO = 10eV(±£) 3 . 

Intriguingly, the dark energy component slightly devi- 
ates from being a pure renormalization constant. Hence, 
as on the upper branch, it has a dynamical character. 
Even more importantly, we find that a component of the 
underlying quantum system, namely, the K^a^ 3 contri- 
bution to the zero-point energy of a state which fulfils 
an approximate KMS condition, shares the redshift be- 
haviour of dark matter. Since this component scales with 



T 3 instead of an inverse volume factor, it does not allow 
for an interpretation in terms of a mass density of non- 
rclativistic weakly interacting particles. So far, the ther- 
mal nature of the quantum state seems to affect only the 
background evolution. It remains to be checked, whether 
its contribution to the energy density is able to repro- 
duce also the clustering properties of dark matter in- 
ferred from astrophysical observations. As a first check, 
we computed (: cj) 2 :)q on a spherically symmetric and 
static spacetime. We estimated the corresponding lowest 
order thermal contribution to the vacuum energy density 
to be proportional to |goo| -3 ^ 2 , where 500 is the coefficient 
of the metric in the static time direction. Since goo can 
be approximated as — (1 + cr 2 ) inside of a spherical body 
|26| . r being the radius, we infer a density profile which 
interpolates between a constant behaviour in the inner- 
most region and a decay with r -3 for large distances from 
the centre. This result stunningly agrees with astrophys- 
ical estimates of dark matter haloes of dwarf spheroidal 
galaxies 27ll28j. 



Outlook. We have applied a quantization scheme which 
is intrinsically suited for curved backgrounds to cosmol- 
ogy. Under the assumption that the underlying quan- 
tum state fulfils an approximate KMS condition at some 
point in the past, we have shown that there exist ho- 
mogeneous and isotropic solutions of the semiclassical 
Einstein equations which are asymptotically stable. The 
exact behaviour of the Hubble function depends on suit- 
able renormalization parameters, intrinsic to the quan- 
tization procedure, which we determined by fitting type 
la supernova data. Interestingly, there are two classes 
of solutions, which are both physically acceptable a pri- 
ori, and which both provide a dynamical interpretation 
of dark energy that reproduces the observed recent ex- 
pansion history of the Universe (at least) as good as the 
ACDM model. 

There are still several open questions to address in fu- 
ture research and here we shall mention only the most 
relevant ones. The first, as already briefly discussed, 
calls us to clarify to what extent the thermal nature of 
quantum states entails an interpretation in terms of a 
new, additional contribution to the dark matter compo- 
nent in the Universe. In order to check whether a full 
quantum description of the complete dark matter sector 
is conceivable, it will be mandatory to study the inho- 
mogeneous fluctuations of the thermal quantum energy 
density and their impact on the formation and growth 
of large scale structures. We also mention the necessity 
first to include interactions and, then, to specify a con- 
crete model. These subsequent steps will bring us closer 
to a more fundamental understanding of the dynamics of 



the Universe. 
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